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1. Introduction 

Prediction of whether a pressurized cylinder will fail catastrophically 
when impacted by a projectile has important applications ranging from 
perforation of airplane's skin by a failed turbine blade to meteorite impact of a 
space station habitation module. This report summarizes the 
accomplishment of one task for a project, whose aim is to simulate 
numerically the outcome of a high velocity impact of pressure vessels. A finite 
element patch covering a vicinity of a growing crack has been constructed to 
estimate the J-integral (crack driving force) during the impact. Explicit 
expressions for the J-integral through the nodal values of displacement, 
strain, and stress have been written. The patch is to be used repeatedly to 
estimate the amount of crack growth during the the time of the impact. The 
resulting crack size is to be compared to an estimated critical crack size for the 
pressurized cylinder. 

A literature search produced a number of papers dealing with evaluation 
of J-integral within finite element environment. Most of the research reports, 
however, present the shape of the finite element mesh only, with no detail on 
node locations. Such information was hard to utilize in the absence of an 
automated mesh generator. As a result, the simplest mesh was chosen for the 
patches, following (2). The same search turned up studies of the accuracy of 
finite element J-integral evaluations as well as the effect of the choice of the 
contour of integration. This provided a rational basis for the choices made in 
the present work. 

A complementary literature search has been done to collect data on 
fracture toughness of 2219 aluminum alloys, since this material property 
enters the employed crack growth criterion. 

The third literature search concerned reports on high- and hypervelocity 
impact studies (both experimental and theoretical) to form a basis for 
comparison with the numerical simulations produced by the entire project. 

Complete computational details and the three literature reviews have 
been left with Rene Ortega. 

2. Circumferential and Axial Patches 

Both patches have the shape of a rectangle with an edge crack mapped 
onto a portion of the cylinder's surface as shown in Fig 1. The finite element 
mesh consists of 8-node isoparametric elements (1). Of these only the four 
which surround the crack tip are distorted, namely, the five nodes neighboring 
the crack tip are placed at the quarter distance from the tip instead of being 
half distance away (see Fig 2). Formulas shown in Fig 1 permit to find the 3D 
coordinates of any node. 

3. J-integral expressions 

J-integral is the following contour integral: 
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crack and ends on the other (see Fig 3). The integral has the meaning of the 
potential energy release per unit crack advance (known as 'the energy release 
rate', or 'the crack driving force’). 

Explicit expressions for the J-integral through the nodal values of 
displacement, strain , and stress have been written for the two contours shown 
in Fig 1. The structure of those expressions is exemplified below for the inner 
contour. 


Eq [1] is rewritten as 

where 

Ii = ^ w dx 2 

r 

and 


[ 2 ] 

[3] 


I 2 = ^ Ojj (dui/Sxi) nj ds M 

r 

The contour is split into five paths r x r 5 (see Fig 1), and the integrals [3,4] 

become the sums of the integrals over these paths: 

Ik = Iki + ••• + Ik5 » k = ® 


As ex am ples, the expressions for I^i and I 2 j through he nodal values of u 2 , £y, 
and Oy are shown here: 

I n - - (h/6) (w 229 + 4W 244 + 2w 255 + 4w 270 + w 281 ) [6] 

I 21 = - (h/6) (f 229 + 4f 244 + 2f 255 + 4f 270 + f 281 ) [71 

where h is the mesh size, the upper indecies refer to node numbers, 


w = o^£y-/2 

f = a n e n + a 12 (du^xj) 

and, as a matter of example, the expression for du^dxj through the nodal 
values of u 2 is shown: 


Ou2/dXi) 244 = (l/2h)(U2 257 + U, 231 - U2 227 - U2 253 ) 

+ (l/h)(U 2 254 + U 2 243 + U 2 228 - u 2 230 - U 2 245 - u 2 256 ) [10] 
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3. Testing of the patch 


To verify the numerical procedures, comparison has been proposed with 
an existing solution for a rectangular plate with an edge crack parallel to the 
clamped edges (4). 

4. Discussion 

The energy release rate and its J-integral representation employed in 
this study corresponded to static (or slowly growing) crack, whereas the crack 
under consideration is a fast growing one. However, it is known that the 
energy release rate for a moving crack is related to the static one as G d y n = 
g(v)G stat , where g(v) is a monotonically decreasing function of the crack 
velocity v which goes from latv = 0to0atv reaching the Rayleigh wave speed 
(3). Therefore employing G stat overestimates the crack driving force and thus 
is conservative when a possibility of a catastrophic failure of the cylinder is 
considered. 

If, nevertheless, the estimates will result in unrealistically large crack 
sizes at the end of the duration of the impact, expressions for dynamic J- 
integrals and their evaluation in finite element environment are available (see 
the literature review). 

Finite element models of elastic-plastic crack growth in the presence of 
both small and large scale yielding are also available in the literature (see the 
literature review). 
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